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ABSTRACT 

The 'Bloated Stars Scenario' proposes that AGN broad line emission originates in the 
winds or envelopes of bloated stars (BS). Alexander and Netzer (1994) established that 
~ 5xl0 4 BSs with dense, decelerating winds can reproduce the observed emission line 
spectrum and avoid rapid collisional destruction. Here, we use the observed properties 
of AGN line profiles to further constrain the model parameters. In the BS model, the 
origin of the broad profiles is the stellar velocity field in the vicinity of the central black 
hole. We use a detailed photoionization code and a model of the stellar distribution 
function to calculate the BS emission line profiles and compare them to a large sample 
of AGN CivA1549, Cm] A1909 and MgnA2798 profiles. We find that the BSs can 
reproduce the general shape and width of typical AGN profiles as well as the line 
ratios if (i) The ionizing luminosity to black hole mass ratio is low enough, (ii) The 
broad line region size is limited by some cutoff mechanism, (iii) The fraction of the BSs 
in the stellar population falls off roughly as r~ 2 . (iv) The wind density and/or velocity 
are correlated with the black hole mass and ionizing luminosity. Under these conditions 
the strong tidal forces near the black hole play an important role in determining the line 
emission properties of the BSs. Some discrepancies remain: the calculated BS profiles 
tend to have weaker wings than the observed ones, and the differences between the 
profiles of different lines are somewhat smaller than those observed. 

Key words: galaxies:active - quasars: emission lines - stars:giant - stars:kinematics 



1 INTRODUCTION 

The observed properties of active galactic nuclei (AGN) lead 
to the conclusion that the broad line emission originates in 
numerous small, cold and dense gas concentrations. These 
objects are labeled "clouds", although their true nature re- 
mains unknown. A basic challenge for any broad line region 
(BLR) model is to specify the physical mechanism that pro- 
tects the clouds from rapid disintegration in the AGN's ex- 
treme environment, or else specify a source for their contin- 



ued replenishment. The bloated stars (BSs) model ( Edward; 



198C|; Mathews 198$ |5coville fc Norman 198g [Penston 1988} 



Kazanas 1989) proposes that the lines are emitted from the 



winds or mass loss envelopes of giant stars. The star provides 
both the gravitational confinement and the mass reservoir 
for replacing the gas that evaporates from the envelope to 
the interstellar medium (ISM). The BS model is further mo- 
tivated by the lack of observational evidence for net radial 
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motion in the BLR (e.g. Maoz et al. 1991, Wanders et al. 
1995). This is consistent with virialized motion in the grav- 
itational potential of the nucleus. 

The BS model was recently studied by Alexander & 
Netzer (1994, paper I), which used a detailed photoioniza- 
tion code to calculate the emission line spectrum of vari- 
ous simple BS wind models. The integrated emission line 
spectrum was calculated by combining the line emissivity of 
the BSs with theoretical models of the stellar distribution 
function (DF). The properties of the models were compared 
to the mean AGN line ratios and to estimates of the BLR 
size and line reddening. The mean observed Lya equivalent 
width was used to determine the number of BSs and their 
fraction in the stellar population and to estimate the col- 
lisional mass loss rate and its effect on the ISM electron 
scattering optical depth. 

The main result of paper I is that there are BS wind 
models that can reproduce the emission line ratios to a fair 
degree with only a small fraction of BS in the stellar pop- 
ulation. However, these models do not resemble winds of 
normal supergiants. The photoionization calculations show 
that the emission line spectrum is dominated by the con- 
ditions at the outer boundary of the line emitting zone of 
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the wind. The successful BS models are those with dense 
envelopes that have small density gradients. In this case the 
wind boundary is set by tidal forces near the black hole and 
by the finite mass of the wind at larger radii. Only ~ 5x 10 4 
such BSs (< 1% of the BLR stellar population) are required 
for reproducing the BLR emission. As a result, the collision 
rate is reasonably small (< IMo/yr). On the other hand, 
lower-density models or those with a steep density gradi- 
ents are ruled out because they emit strong broad forbidden 
lines, which are not observed in AGN. In addition, their low 
line emissivity requires orders of magnitude more BSs in the 
BLR. This leads to rapid collisional destruction, very high 
mass loss rate and very short BS lifetimes. Even if the BSs 
are continuously created, the ISM electron scattering optical 
depth is very large, contrary to what is observed. 

In this second part of the work, we combine our line 
emission calculations with the dynamics of the stellar cluster 
to obtain the predicted line profiles for the BS model. This 
provides a second, powerful test that can, when compounded 
with observed line profiles, put more constraints on this idea. 
Section [] describes the components of the model: the BSs, 
the ionizing continuum, the stellar distribution function and 
the black hole mass. Section ^ describes the line profile cal- 
culations and presents a large sample of observed profiles. 
In section ^ we compare the profiles of various BS models 
to the observed ones. A discussion and a summary of our 
conclusions are presented in section ^. 



2 THE MODEL 

The observed properties of AGN profiles broadly fall into 
the following categories: 

(i) The line width distribution in the AGN population. 

(ii) The profile shapes. 

(iii) The differences between profiles of different lines. 

(iv) Profile asymmetries. 

(v) Profile shifts relative to the systemic redshift. 

A successful BLR model should be able to explain all these 
properties. The BS model, in its current level of detail, does 
not attempt to explain the profile asymmetries nor the rel- 
ative line shifts. Nevertheless, even the remaining profile 
properties set strong constraints on the BS model. We pro- 
ceed on the premise that the asymmetries and relative shifts 
can be eventually explained either within the framework of 
the BS model or by the existence of another line emitting 
component in the BLR. 

The integrated line emission is a spatial average of the 
line emissivity in the BLR volume, weighted by the stellar 
density. Consequently, the line ratios do not impose strin- 
gent constraints on the spatial distribution of the gas. The 
line profile, L u (u is the line of sight projection of the 3D 
velocity v), is a finer probe of the BLR structure, since un- 
like the total emission, it is integrated only over gas which 
has line of sight velocity u. 

The motion of the BSs follows the velocity field of the 
stellar system that surrounds the black hole. The input pa- 
rameters that are required for calculating the line profiles are 
the structure of the BS wind, the black hole mass, the stel- 
lar DF and the distribution of the BSs in the normal stellar 
population. A detailed description of these ingredients can 



be found in paper I. Here we summarize the relevant prop- 
erties and list the new features that were introduced into 
the model. 



2.1 The bloated stars 

There are to date no reliable models for the structure of 
massive winds in the presence of an intense external radia- 
t ion fi eld (see however the re cent w ork by Scoville & Norman 
( 199a| ) and Hartquist et al (1995)). At this stage, our pur- 



pose is to understand the conditions that are required for 
the BSs to be consistent with both the observed BLR prop- 
erties and the properties of the stellar population. We there- 
fore model the BS by a simple structure with no claim to 
hydrodynamical self-consistency and without specifying the 
process that drives the wind. The validity of this approxima- 
tion was demonstrated in paper I by showing that the line 
emission does not depend on the wind internal structure, but 
predominantly on the conditions at the boundary of its line 
emitting zone. The BS is modeled by two components: a gi- 
ant star of radius i?* = 10 13 cm and a mass of M* = 0.8Mq, 
emitting a spherically symmetric wind of radius i? w which 
contains an additional mass of up to M w = 0.2Mq . The wind 
density and velocity are linked by the continuity equation 



Mw = 4nR 2 v(R)N(R) . 



(1) 



where M w = lO _6 M0/yr is the mass loss rate and N is the 
hydrogen number density. The wind structure is parameter- 
ized by a power-law velocity field: 



v(R) =v*(R/R* 



(2) 



where is the velocity at the surface of the star. This family 
of models spans a two-parameter space, (a,v t ), which was 
investigated extensively in paper I. Here we will concentrate 
on the values that best reproduce the observed emission line 
spectrum, a = 1/2 (free falling flow) and v* ~ 10 4 cm/s 
(slow and dense wind, 10 s ^ N <, 10 12 cm -3 ). 

The size of the wind, and hence its boundary density 
is determined by three competing physical processes. Tidal 
disruption of the outer layers of the wind by the black hole 
limits the size of the wind to 



Rvm(r) = X tida i(M*/M bh ) 1/3 r , 



(3) 



where Mbh is the black hole mass, r is the distance from the 
black hole and Xydai is a factor of order unity (we assume 
-Xtidai = 2). We will use the term "tidally- limited" to des- 
ignate such BSs. The finite mass in the wind sets an upper 
limit on the wind size, fl maBS , 

fRn 



dR = M w 



R* 



(4) 



We will use the term "mass-limited" to designate such 



BSs. The third process is Comptonization (Kazanas 1989) 



whereby the central ionizing continuum heats the outer lay- 
ers of the wind, ionizes them completely and reduces their 
optical depth to zero. The radiation penetrates into the 
denser parts of the wind until the density rises above a crit- 
ical value Ac, and an ionization equilibrium is established 
at the Comptonization radius Rc- The emission line spec- 
trum is determined by the ionization parameter, U, which 
is essentially the ratio of the ionizing flux and gas density 
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at the gas surface. The results of paper I show that the gas 
at the Comptonized wind boundary has a high U. This re- 
sults in strong emission of unobserved forbidden lines, and 
therefore this process must be suppressed. The high gas den- 
sity and small density gradient of the slow and dense wind 
models achieve this suppression by having Rc > -Rtidai or 

Rc> R mass - 



2.2 The galactic nucleus 

The stellar distribution function that we use i s base d on the 
numeric results of Murphy, Cohn & Durisen (1991) (hence- 



forth MCD), which follow the evolution of a multi-mass co- 
eval stellar cluster in the presence of a central black hole. 
The black hole mass grows as it accretes mass from the stars, 
both directly, by tidal disruption, and indirectly, from mass 
loss in the course of stellar evolution and of inelastic stellar 
collisions. The MCD calculations assume spherical symme- 
try and an initial Plummer DF with a seed black hole of 
1O 4 M . 

In this work we use MCD model 2B, which has an ini- 



tial central density of 7xl0 7 M Q /pc 3 



3xlO 8 M within 



the inner 1 pc and a total mass of A/ sys = 8.5x10 Mq 
within 100 pc. The luminosity is due to spherical accre- 
tion with rest-mass to luminosity conversion efficiency of 0.1. 
The luminosity reaches a peak of 2.4xl0 46 erg/s at 4xl0 8 
yr. The black hole reaches a mass of 5.5x10 8 Mq after 15 
Gyr. We study this system at two epochs. The first is at 
to = 3xl0 8 yr, when M hh = 8xl0 7 M Q , L ion = 7xl0 45 
erg/s and the luminosity is Eddington-limited. This young 
system, which was the one studied in paper I, has the lu- 
minosity of a bright quasar. The second is at to = 10 9 yr, 
when M bh = 1.9xl0 8 M Q , L ion = 3.6xl0 44 , the luminosity 
is sub-Eddington and determined by the mass loss rate from 
the stellar system. This evolved system has the luminosity 
of a luminous Seyfert 1. 

MCD assume, for simplicity, that the velocity DF can 
be approximated by the Maxwell-Boltzmann DF 



DF(r,v,t) 



3 \ 3 / 2 n*(r,t) 



,2-71"/ vo(r, t)' c 



exp 



vo(r, t) 



(•») 



where n* is the stellar density and vo the r.m.s velocity. 
Thus, the velocity field is fully specified by vo(r, t), which 
MCD calculate numerically by following the time evolution 
of the system. 

Our calculations use an approximate n* (r, which is 
based on the density function of the O.8M0 mass bin stars 
(that being the mass bin closest to IMq in MCD figure 6, 
see paper I). The density function no.s(r) is given in MCD 
for several values of to. We interpolate it to any to and then 
approximate n* by 



n*(r,to) = n . s (r,t ) 



Mbh 



MqAtt J r' 2 no.sdr' 



(6) 



that is, by approximating that the ratio between the to- 
tal stellar density and that of the O.8M0 mass bin stars, 
(~ 4), does not depend on r and that the mean stellar 



3 Note that here, as opposed to paper I, we use to designate 
the total stellar population, BSs included. 



mass is IMq. For vo(r, t) we use a semi-analytic approxi- 
mation, which is derived from the following arguments. We 
model the stellar core at times to = 3x 10 8 and 10 9 yr, which 
are both earlier than the collisional time-scale of the initial 
Plummer distribution, t c — 1.2 Gyr. The relevant dynamical 
length scale in the nucleus is the black hole dynamical ra- 
dius, rdyn = GMbh/v? nit , where vixut is the initial Plummer 
r.m.s velocity (v init = 1360 km/s for MCD model 2B). At 
r <C r c i yn , the initial DF evolves rapidly due to the frequent 
stellar collisions in the deep potential well, so that by to the 
velocity field is dominated by the black hole. Consequently, 
only the black hole mass, rather than the detailed form of 
n*(r, t), is required for calculating no. The spatial redistribu- 
tion of the stellar mass at r <C rd yn , whether by the evolution 
of the DF or by the growth of the black hole, does not af- 
fect the potential at r > rd yn because the mass distribution 
is spherically symmetric. Consequently, if to < t c , the ini- 
tial Plummer DF is almost unchanged except for the small 
evolutionary mass loss. The MCD models assume that the 
initial mass function is independent of r, and therefore the 
evolutionary mass loss is reflected only in a small, spatially 
homogeneous decrease of the initial stellar density. 

At small r , the potential of the black hole dominates and 
v'o = yGMhh/r, where 7 is of order unity. The width of the 
line profiles strongly depends on 7 and it is therefore impor- 
tantto estimate it as reliably as possible. Cohn & Kulsrud 
(jl97§) fi nd th at 7 ~ 1.2 at 0.001 < r/r dyn < 0.1. Bahcall 
& Wolf (Ejffll similarly find that 7 J> 12/11 at r < r dy „ 



and that it increa ses by up to a factor of 2 as r decreases. 
David et al (1987) and MCD (B. Murphy, private communi- 



cation) find that 7 approaches 2 as r decreases. This result 
can be explained qualitatively in terms of the orbits of the 
stars near the black hole. If these stars were in bound cir- 
cular Keplerian orbits, 7 would have been 1. However, such 
bound stars are quickly destroyed by stellar collisions and 
tidal disruption. The only stars that can exist at such small 
radii are those that spend there only a very short fraction 
of their period and have high enough angular momenta to 
avoid falling into the black hole. Such orbits are marginally 
bound, nearly parabolic Keplerian orbits, which have 7 = 2. 
At larger radii the stellar destruction proceeds more slowly, 
the fraction of bound stars increases and therefore 7 de- 
creases. 

The r.m.s. velocity of the initial Plummer distribution 



wo = 



*(r) 



where ^ is the negative of the gravitational potential 
|jy _ GM SYS 



(7) 



(8) 



and ro = 1 pc is the core radius (Binney & Tremaine 1987). 
We propose to approximate the r.m.s velocity by the sim- 
plest expression that approaches the correct limits at small 
and large r, 



2, > G(M. 
v (r) 



+ l{r)— — 



(9) 



where y(r) is a smooth function that approaches 2 at small 
radii and 1 at large radii. 

We approximate 7(7") by assuming that dynamically, 
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the stellar population near the black hole is composed of 
two types of stars. A fraction /b of the stars are bound to 
the black hole in circular orbits and have 7t, = 1, while the 
remaining f u = 1 — /b of the stars are unbound, in parabolic 
orbits with y u — 2. We further assume that the velocity dis- 
tributions of the bound and unbound stars are Maxwellian 
with Wq b and vq u of the same form as equation ^, with 7t, 
and 7«, respectively. It then follows that the total popula- 
tion has Vq — fbv'o b + /««o,u and that 



j(r) = ffyb + ful-u 



h 



(10) 



We assume that initially all the stars are bound. With time, 
inelastic collisions and tidal disruption by the black hole de- 
plete the fraction of bound stars. The MCD results show that 
inelastic stellar collisions is the dominant mass loss process 
for the high density stellar cores which are studied here. We 
therefore neglect the tidal disruptions as a mass loss process. 

fbito) is equal to the probability of a bound star to 
survive up to time to- At each inelastic encounter, which 
occur at a rate q c , the star loses on average a fraction k 
of its mass. By approximating that the fractional collisional 
mass loss rate, nq c , is constant in time and equal to its value 
at to, we obtain 



f b = exp(-Kq c t ) . 



(11) 



The collisional rat e for equal mass and radiu s stars with a 
Maxwellian DF is ([Binney fc Tremaine 1987[) 



Qc(to) = 16y -n+v Rl 



1 + 



(12) 



where R+ is the stellar radius, v osc is the escape velocity from 
the star and the second term in the brackets is the Safronov 
number, which takes into account the enhancement of the 
collision rate by gravitational focusing. For k we use the 
approximate form suggested by MCD, 



c/«0 



= 0.1. 



(13) 



Both q c and k are functions of 7 (through vo) and therefore 
7(7-) (equation |lo| ) is the root of the equation 



7 = 2- exp(-K(7)g c (7)t ) , 



(14) 



which is solved for 7 numerically. 

Figure ^ compares the approximated vo to that calcu- 
lated by MCD for model 2B^] at two epochs (B. Murphy, 
private communication) . The approximation tends to under- 
estimate the velocity at larger r and slightly over-estimate 
it at small r. These trends appear also for the younger or 
less massive MCD models. However, in all these cases the 
fractional error in the approximation is less than 10%, which 
we consider to be well within the joint uncertainties of the 
BS model and the MCD results. 

Finally, we note that BSs are expected to undergo more 
collisions than the smaller main sequence (MS) stars. It 
is therefore possible that the BSs will not share the same 
velocity field as the general stellar population but are all 
marginally bound with 7 = 2. Figure [l] shows that this may 

^ The numeric results displayed in Fig. [l] were calculated for an 
initial mass function in the range 0.12 to 12Mq instead of the 
original range of 0.3 to 3OM0 in MCD. 
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Figure 1. The analytic approximation of the r.m.s stellar veloc- 
ity, vo, compared to the MCD results. Solid line: The approximate 
vq (equation ^) of the main sequence stars. Dashed line: The ap- 
proximate vo of the marginally bound (7 = 2) stars only. Open 
circles: vo calculated by MCD for model 2B at to = 6.3 xlO 7 yr, 
when M b h = 6.0xl0 7 M Q . Full circles: v calculated by MCD for 
model 2B at t = 5.3x10 s yr, when A/ bh = 1.3xlO 8 A/ . (Cour- 
tesy of B. Murphy) 



result in an increase of up to ~ 20% in the BS r.m.s veloc- 
ity in young systems when to is much smaller than the MS 
collisional time-scale, but less so in more evolved systems. 

2.3 New features of the BS model 

2. 3. 1 The fraction of BSs 

In paper I, we assumed that the fraction of the BS in the 
stellar population, /bs, is constant and independent of r. 
This is obviously only a crude approximation. BSs are not 
observed in normal stellar environments, such as are, pre- 
sumably, the host galaxies of AGN. This implies that /bs 
must approach zero at large r. The results of section ^ be- 
low indicate that the form /bs oc r~ 2 may be necessary for 
reproducing the shape of the observed profiles. We therefore 
model here also the case of 

/bsM = min(/ max ,/ max (jy) ) , (15) 

where / max < 1 is an upper bound on the fraction of BSs 
in the stellar population and r/ is the distance where the 
BS fraction reaches this bound. / max is a free parameter 
of the BS model. We will make the arbitrary assumption 
that / max ~ 1/4, the fraction of the 0.8Mq mass bin stars 
in the stellar population. For the model discussed below, 
vo(rf) > 10000 km s _1 . Since BSs with such a high velocity 
contribute very little to the line emission, for all practical 
purposes /bs <x r~ 2 . 



2.3.2 Collisional mass loss 

Collisional mass loss sets a strong constraint on BS wind 
models. Reliable estimates of the mass loss rate are espe- 
cially important for BS distributions that are concentrated 
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towards the black hole. Here we refine the simple estimates 
used in paper I for the collisional mass loss rate per volume 
by using equations ^| and ^ and by taking into account the 
exact amount of mass in the wind, rather than the typical 
M w = 0.2M©. Thus, for BS-BS collisions 



PBScoiiM = /BsM n *( r )[<7c,*re*M* + 
(<?c,w — q c ,*)K w M v/ (R w )] , 



(16) 



where the collisions of the giant star and of the wind are 
treated separately, q c ^ and re* are calculated with M+, R* 
and q c ,w and calculated with M* + M„(R„), i? w . 

Similarly, for the MS-MS collisions 

PMScoiiM = [1 - /BsW]n*M?cKM , (17) 
where q c and re are calculated with 1 Mq and 1 Rq. 



drr Li(r) 



dv\v\DF(r,v) 



(19) 



where the 5 function is translated into the integration range 
{ti,T2), which is the sub-interval of (rj n ,r out ) where \u\ < 
Vmax(r). For the assumed Maxwell-Boltzmann velocity field 
(equation |^) v msx (f) — * oo and the line profile is 



L u — 



drr 2 LAr) *\ } exp 
v (r) 



Vo(r) 



(20) 



Equation |20|, together with the Li (r) results of the photoion- 
ization calculation and the approximate Vo(r) (equation ^), 
are used below to predict the line profiles of several different 
BS models. 



2.3.3 Reddening correction 

The amount of intri nsic reddening in AGN spectr a is still 
an unresolved issue ( MacAlpine 1985 ; Netzer 1990 ). We as- 
sume the "screen hypothesis", whereby the dust lies just 
outside the line emitting region and intercepts all the ob- 
served spectrum, continuum as well as lines. In paper I we 
estimated the extinction coefficient, Eb-v, from the differ- 
ence between the calculated and the mean observed Lya/H/3 
ratio. We then corrected all the other lines accordingly. This 
procedure does not allow for the fact that there is actually 
a range of observed Lya/H/3 ratios. It is also vulnerable to 
possible inaccuracies in the Lya and H/3 radiative transfer 
method that is used by the photoionization code. This was 
not a serious problem in paper I, since the models that suc- 
ceeded in reproducing the emission line spectrum required 
only a small reddening correction, while the other models 
were ruled out by excessive collisional mass- loss. Here we 
relax the reddening criterion and determine -Eb-v within 
the range to _EB-v(Lya/H,3) by a % 2 fit of all the calcu- 
lated line ratios and their mean observed values (table 1 in 
paper I). 



3 THEORETICAL AND OBSERVED LINE 
PROFILES 

3.1 Theoretical line profiles 

The basic output of the photoionization calculations is 
Lt(r), the luminosity in a given line from a single BS at 
a distance r. The numerical procedure for calculating Le(r) 
for a given BS wind density structure, chemical composition 
and the ionizing flux spectrum is as described in paper I. 

L u , the differential line luminosity per line of sight ve- 
locity u, is given by 

L u = / drr 2 / dQ, r I dvv 2 I dfl v x 

Jri n Jo Jo Jo 

L t (r)BF{r,v)5(vfir-u), (18) 

where DF(r,v) is the BS distribution function and D ma!! (r) 
is the maximal stellar velocity at distance r from the black 
hole. For an isotropic stellar distribution with an isotropic 
velocity field, DF(r,«) = DF(r, v), L u reduces to 



3.2 Observed line profiles 

The calculated line profiles are compared to a sample 
of 87 C ill] A1909 profiles, 140 C IV A1549 profiles and 72 
MgllA2798 profiles, for which we measured the full line 
widths (FW) at 1/4, 1/2 and 3/4 of the maximum height. 
This sample was compiled from several sources. Ill of 
the CivA1549 lines are fro m an absorption line survey of 



123 high luminosity AGN ( jargent, Boksenberg fc Stei- 



tlcl 198 ISargent, Steidel fe Bokscnbcrg 198S| ; |Sargcnt~ 



Steidel 1995), which were selected by apparent magnitude 
and redshift. The QSO lie primarily in the redshift range 
z ~ 1.9 to 3.5 and the in the absolute magnitude range 
Mv ~ —32 to —27. The continuum subtracted spectra of 
thi s sam ple, which were measured here, are from Wills et 
al 



(199: 



The Cm]A1909 and MgllA2798 and 29 of the 
CivA1549 profiles are from a sample of 42 radio-selected 
QSO and 50 radio-quiet QSO , which were selected b y UV 



excess or slitless spectroscopy ( pteidel fc Sargent 1991 ). The 
QSO lie primarily in the redshift range z ~ 0.9 to 2.2 and 
in the rest-frame monochromatic absolute magnitude range 
M(A2200A) r 30 to -26. The continuum subtracted spec- 
tra of this sam ple, wh ich were measured here, are from 
Brotherton et al (1994a). All the continuum subtracted spec- 
tra were kindly supplied by M. Brotherton. 

We are interested in studying the shape of the line pro- 
files, which requires information on the line width at several 
heights. A common procedure for measuring the FW of a 
spectral line is to fit the profile to a superposition of sev- 
eral Gaussians and measure the width of the smooth fitted 
function. This procedure is also used to deblend adjacent 
spectral lines from the profile. Since we are interested in the 
large sample statistics rather than the properties of a spe- 
cific AGN, and since the art of multi-Gaussian fitting has 
a strong subjective element in it even under optimal con- 
ditions, we decided to adopt a simpler procedure. We rely 
on the fact that all lines of interest are much stronger than 
neighbouring lines and do not attempt any deblending or 
fitting, but simply located "by eye" the profile peak and un- 
derlying continuum level. In cases where the blended line was 
clearly visible, we adjusted the FW "by eye" accordingly. In 
cases where the spectrum was very noisy, we smoothed it by 
averaging the flux in overlapping windows of n consecutive 
pixels (usually n = 5). 

The measurement of the MgllA2798 FW is compli- 
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Figure 2. The average observed CivA1549, Cin]A1909 and 
MgllA2798 profiles, weighted by the integrated flux. The 
Mg II A2798 profile is shown both with and without the assumed 
Fell contribution. 



cated by the fact that this line is blended with a number 
of Fe II lines, whic h add to the Mg ii A2798 pro file wide and 
flat broad wings flWills, Netzer fe Wills 1985|). Since it is 



not our intention to do here a careful deblending of the 
Mg II A2798 profile, we measured the Mg II A2798 twice, once 
with these wings and once without them. In the latter case, 
only 68 profiles were measured since the others lacked dis- 
tinct wings. These two types of measurements, which we 
designate by Mg II A2798+Fell and Mg II A2798-Fe II , re- 
spectively, represent the two extreme assumptions that ei- 
ther the wings are all Mg II A2798 emission or that only 
the line core is Mg II A2798 emission. Thus, a BLR model 
that predicts FW(Mg n A2798) < FW(Mgn A2798-Fe n ) or 
FW(MgnA2798) > FW(Mgn A2798+Fen) is inconsistent 
with the observations. 

We compared our FWHM resu lts with the Gaus- 
sian fitresults of Wills et al (1993) and Brotherton et 
al ( 1994a[ ) by calculating the mean fractional difference, 
(FWHM by eye - FWHMg auS s)/FWHMgauss, and its rms 
scatter. For the CivA1549, Cm]A1909, Mg II A2798+Fe II 
and Mg II A2798— Fe II profiles the mean fractional differ- 
ences are -0.009T0.06, 0.05T0T5, -0.01T0T and -0.26T 
0.13, respectively. The correspondence of the two meth- 
ods for the CivA1549 and Mg II A2798+Fe II lines is very 
good. The relatively larger over-estimation of the C ill] A1909 
FWHM is due to the blended AlmA1859 and the big dif- 
ference in Mgll A2798— Fen is due to the fact that Wills et 



al ( |1993j ) and Brotherton et al ( |l994a| ) did not deblend the 
Fe II contribution to the Mg II A2798 profile. 

Figure ^| shows the sample averaged profiles, where the 
profiles are weighted by their inverse integrated flux. It 
should be noted that other weighing methods (e.g. by the 
peak flux), or other averaging procedures (e.g. averaging the 
blue and red half widths at x% of the maximum, for x = 
to 100%) yield very different results. This ambiguity indi- 
cates that the average profile is not a very useful quantity 
and that the calculated profiles should be evaluated by the 
cumulative distribution function or on an object to object 



basis. Figure |§] shows the FW cumulative distribution func- 
tions of the three lines at 1/4, 1/2 and 3/4 maximum. 



4 RESULTS 

We investigate the line profiles of the BS model by studying 
a succession of three BS wind / stellar DF models, whose 
properties are summarised in table [jj. Model A is similar 
to that studied in paper I. We show below that in order to 
reproduce the observed line profiles and ratios, some changes 
in the model assumptions are required. These changes are 
introduced in two steps, in models B and C. 



4.1 Model A 



We begin by considering the successful BS model of paper 
I (hereafter model A, see table ^), which has the BS wind 
parameters a — 1/2 and = 5xl0 3 cm/s, an AGN model 
of M bh = 8x 10 7 Af Q and L ion = 7x 10 45 erg/s at t = 3x 10 8 
yr and /bs = const. A free parameter of the model is the 
BLR size, r out . As was shown in paper I, the line spectrum 
depends only weakly on r out in the range 1/3 pc to ~ 50 
pc. This is not the case for the line width, which decreases 
significantly with r ou t- We will assume that the BLR outer 
limit scales as r ou t ~ (Lion/10 46 erg/s) 1/2 pc. This is mo- 
tivated both by the line reverberation results that suggest 
that the average BLR radius is ~ 0.1(L io n/10 46 erg/s) 1/2 pc 
(Netzer 1990) and by the possibility that the line emis- 
sion is suppressed by dust at th e dust sublimation r adius 
r sub ~ 0.3(L io n/10 46 erg/s) 1 / 2 pc ( |Netzer fe Laor 1993[ ). We 
therefore set r ou 



1 pc for model A. 
Figure ^ shows some profiles calculated for a BLR size 
r ut = 1 pc. Three properties are immediately apparent: The 
profiles are very narrow (FWHM ~ 2200 km/s), the wings 
are very weak and the different lines have almost identi- 
cal profiles, unlike the observed ones (Figure ^). Table || 
shows that only a few percent of the objects in the sample 
have FWHM as narrow as those calculated, and none have 
such narrow FW1/4M. We explored also the possibility of 
reducing r out to 1/3 pc. This increased the CivA1549 FW 
to 1877 (> 29% of the sample), 2982 (> 14%) and 4422 
(> 5%) km/s for the FW3/4M, FW1/2M and FW1/4M, re- 
spectively. However, this modest increase in FW1/4M comes 
at the price of an excessive BS-BS collisional mass loss rate 
of 4.5Mq /yr. As will be discussed below, the emissivity of 
tidally-limited BSs is approximately independent of r. De- 
creasing the BLR size therefore means that the same num- 
ber of BS must share a smaller volume and this leads to an 
increased collision rate. 



4.2 Model B 

A possible remedy for the width problem is to model a more 
evolved AGN with a more massive black hole. The MCD 
models indicate, that apart for a very short initial epoch of 
Eddington limited accretion, the luminosity decreases with 
time as the available mass is exhausted (MCD figure 10). 
The resulting decrease in r out also helps to broaden the 
profiles. We therefore consider an AGN model at to = 1 
Gyr, when M hh = 1.9xl0 8 M Q , L ion = 3.6xl0 44 erg/s and 
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Figure 3. The cumulative distribution functions of the observed FW at 1/4, 1/2 and 3/4 maximum for CivA1549, Cm] A1909 and 
MgnA2798. The QSO number density per FW interval is proportional to the density of vertical lines. The Mgn A2798 cumulative FW 
distribution is based on the Mg II A2798+Fe II measurements. The Mgn A2798— Fen distribution is outlined by the bold line. 



Tout = 1/4 pc (hereafter model B, see table As is ex- 
plained in the appendix, changing Af D h and Li on from the 
values of model A to those of model B worsens the fit of the 
calculated line ratios to the observed ones, unless accompa- 
nied by a suitable modification of the wind parameters. In 
paper I, we systematically searched the BS wind parameter 
space for the values that best reproduce the observed line 
ratios. Here, we use instead the scaling prescription given in 
the appendix, which shows that the line spectrum of model 
A can be approximately reproduced by modifying the wind 
velocity to = 8xl0 4 cm/s. For the continuum spectrum 
we use the same spectral distribution as used for model A, 
scaled down to the lower Li on of model B. 

The calculated BS-BS collisional mass loss rate of 
model B is only OAMq/yt, in spite of the small BLR size. 
This reflects the fact that most of the mass loss is from col- 
lisions of the winds, and that the BS wind density is propor- 
tional to M-w/v* (equation^)). Since J? t idai is a weak function 
of Mhh) the BS wind mass in model B is ~ 16 times smaller 



than that of model A, and so are the respective mass loss 
rates. Figure ^] and table |^ show that the profiles of model 
B are indeed wider than those of model A, but the profiles 
of the three lines are almost identical and the wings are too 
weak. 

So far, we assumed that the fraction of BSs in the cen- 
tral cluster, /bs, is independent of r. There are, however, 
good reasons to modify this approach, both from theoreti- 
cal considerations and from careful examination of t he lin e 
wings. As noted by Penston, Croft, Basi & Fuller (1990), 



the far wings of many observed profiles are well fitted by an 
inverse quadratic form of Fx oc 1/(A — Ao) 2 and show that 
this is consistent with line emission from an ensemble of 
clouds of constant cross section, moving on parabolic orbits. 
Here, we extend their analysis to the BSs and show that the 
discrepancy between the observed wing shapes and those of 
model B can be explained in terms of the line emissivity 
from marginally bound BSs. 
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Model A 


Model B 


Model C 




5x 10 3 cm s _1 


8xl0 4 cm s" 1 


8xl0 4 cm s" 1 




8.0xl0 7 M Q 


1. 9x10 s M 


1.9xlO 8 M 


to 


3xl0 8 yr 


10 9 yr 


10 9 yr 


-^ion 


7xl0 45 erg s 1 


3.6xl0 44 erg s -1 


3.6 Xl0 44 erg s -1 


nil 


0.001 pc 


0.001 pc 


0.001 pc 


r ou t 


1 pc 


0.25 pc 


0.25 pc 


No. of BSs 


1.5xl0 5 


4.3xl0 4 


4.7xl0 4 


Fraction of BSs 


4.7xl0 -4 


2.4xl0 -3 


2.6xl0 -3 


/bs 


oc r° 


oc r° 


OC T z 


M evo \ 


0.15M Q yr" 1 


4.3xlO~ 2 M yr" 1 


4.8xlO- 2 M Q yr" 1 


A^BS coll 


3.2M Q yr" 1 


0.37M Q yr" 1 


O.59M yr" 1 


A^MS coll 


2.8xlO" 2 M Q yr" 1 


I.lxl0- 2 M Q yr" 1 


l.lxlO~ 2 M yr" 1 


Electron scattering f 


< 3.8xl0~ 4 yr" 1 


< 1.8X10" 4 yr" 1 


< 0.08 yr" 1 


Eb-v 


0.066 


0.000 


0.096 


hya flux weighted radius 


493 Id 


209 Id 


25 Id 



Table 1. The parameters and calculated properties of models A, B and C. The growth rate of the electron scattering optical depth, f, 
is estimated under the assumption that all the mass accumulates in the BLR and is neither accreted nor ejected (see paper I). 



Model A Model B Model C 



FW 


CivA1549 


C in] A1909 


Mgn A2798 


CivA1549 


Cm] A1909 


MgoA2798 


CivA1549 


C in] A1909 


Mgn A2798 


3/4M 


1337 


1407 


1252 


2398 


2426 


2366 


2854 


2958 


2729 


% 


(4) 


(8) 


(1/5) 


(49) 


(38) 


(54/91) 


(61) 


(47) 


(76/96) 


1/2M 


2090 


2202 


1952 


3734 


3789 


3709 


4561 


4782 


4319 


% 


(2) 


(2) 


(0/3) 


(33) 


(15) 


(41/84) 


(48) 


(31) 


(62/89) 


1/4M 


2988 


3164 


2783 


5372 


5464 


5813 


6928 


7344 


6424 


% 


(0) 


(0) 


(0/0) 


(12) 


(2) 


(16/86) 


(32) 


(23) 


(20/90) 



Table 2. Comparison of the calculated FW to the observed cumulative FW distribution. The FW are in km/s. The numbers in brackets 
below each FW entry are the percentage of objects in the sample with FW less or equal to the calculated one. For the MgllA2798, 
the percentage on the left is based on the Mg II A2798+Fe II FW measurements, and that on the right on the Mg II A2798— Fe II FW 
measurements. 

very nearly constant .This, however is not true for the in- 
dividual lines (Fig. ||), whose emissivity depends on the gas 
density, ionization level and optical depth. Nevertheless, let 
us make the approximation that Le ~ const, (a = 0) and as- 
sume also that /bs = const. (6 = 0). Marginally bound stars 
have n* (r) oc r" 1/2 (c = -1/2) (e.g. Penston et al 1990). 
The profile wing is therefore expected to fall off sharply as 
L u oc tt -6 . Figure [| shows that the calculated profiles of 
model B are well fitted by this expression at high it, thereby 
also justifying the approximation of the DF by the 8 func- 
tion. 



4.3 Model C 

4-3.1 Line profiles 

Equation ^ and the observed profiles suggest that in or- 
der to reproduce the observed inverse quadratic wings, it is 
necessary to modify some of the model assumptions, so that 
a + b + c = 5/2. Changing the BS emissivity by modifying 
the wind structure will change the emission line spectrum, 
as shown in paper I. The r dependence of the stellar DF is 
based on a well established result. The only remaining free 



Let us assume for simplicity that the Maxwellian DF 
(equation |H|) can be approximated by 

^(r,») = fa(rW(r) f ' B ^ , (21) 
where v 2 , — 2GMbh/r. It then follows from equation |l9] that 



L u = 2n *r 2 Li(r)yBs(rK(r)t;o 1 (r), (22) 



where the upper integration limit is defined by r-z = 
min(2GA/bh/^ 2 , fout) and 3> r\. As explained earlier, 
marginally bound stars have t>o( r ) oc r -1 / 2 . Assume that 
Li{r) oc r~ a , /bsM oc r~ and n*(r) oc r~ c . It then follows 
that 

L u oc r y 2 -( a+b+ ^ . (23) 

At the profile wings u is large enough so that T2 < r ou t, 
rau" 2 and the profile shape is 

L u OC ^+2(a + 6+ c) . (24) 

The tidally-limited BSs are dense and optically thick 
enough, so that their total emissivity is roughly propor- 
tional to the product of their geometrical cross-section, 
7ri? 2 dal oc r 2 (equation []) and the central irradiating flux, 
F oc r~ 2 , and is therefore constant. Figure ^| shows the 
BS line and total emissivity of model B. The total emit- 
ted heat (including all lines and diffuse continua) is indeed 
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Figure 4. The theoretical C IV A1549, C III] A1909 and Mg II A2798 
profiles for models A, B and C. The model B numerical profiles 
arc fitted to the analytical approximation L u oc ti — 6 (bold dash- 
dotted line). The model C numerical profiles are fitted to the 
analytical approximation L u oc u~ 2 (bold dash-dotted line). 



Figure 5. The line emissivity of a single BS, L(, for models B 
and C in various lines, as function of the distance from the black 
hole. The total heat includes all the reprocessed radiation emitted 
in lines and continuum radiation. 



parameter is /bs- If we assume /as oc r then the profile 
wings will have the required shape^. 

Model C (table ^) has the same parameters as model B 
but with /bs oc r~ 2 (equation^). Figure^] and table ^show 
that with this choice, the calculated FWHM are represen- 
tative of the observed sample. However, the percentiles de- 
crease from the FW3/4M to FW1/2M to FW1/4M, indicat- 
ing that the calculated wings are still somewhat weaker than 
the observed ones. The calculated MgnA2798 FW3/4M is 
now in fact wider than more than 76% percent of the sample 
profiles. The fit to L u oc u~ 2 , shown in figure^ is not perfect. 
This is due to the fact that Lt of those lines is not constant. 
We have examined the fit to the calculated Lya profile and 
found that for this line, whose emissivity is closer to being 
constant with r, The fit of the far wings to L u oc u~ 2 is very 
good. 

The FWHM ratios of the calculated lines of model C are 
Cm] A1909/CivA1549 = 1.05, Mgn A2798/Civ A1549 = 
0.95 and Mg II A2798/C ill] A1909 = 0.90. These values can 
be compared to the observed values quoted in Brotherton et 
al (1994), 1.21 ± 0.04, 0.97 ± 0.06 and 0.85 ± 0.05, respec- 
tively. It is apparent that the model C in] A1909 tends to be 
narrower than the observed one, as is also indicated by the 
percentages in table |^. Figure ^ demonstrates that there are 
AGN which have profiles that are very similar to those of 
model C. 



4-3.2 Line ratios 

Having established a good agreement between the observed 
and calculated profiles, we now further test model C by re- 
turning to the issues studied in paper I, namely the line 
ratios and the required number of BSs. The integrated line 



II /bs oc r 2 corresponds to a differential geometric covering fac- 
tor that falls off as oc r~ 1//2 , since i?tidal oc r and n* oc r~ 1 / 2 . 
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Figure 6. Examples of observed profiles (solid line) which are well fitted by model C (dotted line). 
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Figure 7. The unreddened (dashed line) and reddened, Eb—v = 
0.096, (solid line) spectrum of model C compared with an average 
observed QSO spectrum (dotted line). 



spectra of models A and C are compared to the mean ob- 
served one (paper I, table 1) in figure [| (Note the similarity 
between the line ratios of models A and C, which demon- 
strates that the scaling prescription of the appendix works). 
The MgnA2798 and NvA1240 deficiencies, which are dis- 
cussed in paper I, remain a problem. Model C has a slightly 
high [O ill] A4363 emission, which is due to the fact that the 
scaling only deals with the tidally-limited BSs. At large r the 
BS are mass-limited. The lower wind gas density of model 
C implies a larger i? maS s and consequently a lower density, a 
larger [/(i? maS s) and a stronger emission of forbidden lines. 
The unreddened and reddened spectra of model C are com- 
pared to a mean QSO spectrum in figure |^. 

Model C is quite successful in reproducing the BLR 
emission properties. The line ratios, line profiles and BLR 
size, as estimated by the flux weighted radius of the Lya. 



[FeXI] A7892 
[NeV] A3426 
[Oill] A5007 
[On] A4363 
CIII A977 
Hell A 1640 
Hel A5876 
Mgll A2798 
CIII] A1909 
OVI A 1035 
NV A1240 
CIV A 1549 



0.2 0.4 0.6 0.8 1 12 1.4 1.6 1.8 2 
L(line) / L(Lya) 

Figure 8. The integrated line ratios of models A (A) and C (□), 
expressed relative to the Lye? luminosity and compared to the 
observed range (paper I, table 1). The mean observed ratios are 
indicated by the black dot and the ranges by the error bars. 
E B -v(A) = 0.066 and r out (A) = 1 pc. E B -v(C) = 0.096 and 
rout(C) = 1/4 pc. 



The total mass loss rate is consistent with the theoretical 
estimates of the accretion rate and the resulting electron 
scattering optical depth growth is consistent with the obser- 
vations of rapid X-ray variability. 



5 DISCUSSION AND CONCLUSIONS 

The results of paper I indicate that in order to reproduce 
the integrated emission line spectrum, the BSs should have 
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dense, decelerating winds, whose boundary is not fixed by 
comptonization. The BLR line profiles are another observed 
property that the BS model has to reproduce. We find that 
this imposes additional constraints on the model. 

The FW statistics of the observed profiles indicate that 
the high Li on /Mbh models that were successful in reproduc- 
ing the emission line spectrum have too narrow FWs. Ob- 
served AGN with such narrow broad lines constitute only a 
marginal percentage of the sample. Another potential prob- 
lem with these models is the rather high collisional mass 
loss rate. The narrow calculated FWs indicate that either 
Mbh should be larger and/or the BLR size smaller. The 
MCD models of the AGN evolution link the increase in 
the black hole mass with a decrease in the luminosity (in 
the post-Eddington limited accretion epoch). The change in 
the luminosity affects the model properties in two different 
ways. First, it may decrease the BLR size, since it is natu- 
ral, although not necessary, to assume that the luminosity 
determines the BLR size by creating a dust-free inner re- 
gion which can emit lines. Second, the reduced BLR size 
and increased black hole mass make the tidally-limited BSs 
the main contributors to the line emission. We find that, 
in order to obtain the observed line spectrum from tidally- 
limited BSs, it is necessary to scale M w /«* in proportion 
to Li on /V Mbh (see the appendix). This in turn reduces the 
wind mass and allows the BSs to survive the increased col- 
lisional rate in the smaller, faster moving BLR. The lower 
Lion/Aft* models have much wider profiles, which are con- 
sistent with the observations. However the far wings are very 
weak and the profiles of the different species are very similar 
to each other, contrary to the observed situation. 

Although the profile widths of model A are much 
smaller than those of the sample presented here, we note 
that there exists a sub-population of narrow line Seyfert 1 
galaxies, which are characterized by soft X-ray spectra and 
comprise about 10% of optically selected Seyfert Is and 15- 
50% of soft X-rays selected Seyfert Is (Boiler, Brandt & 



Fink IE 95 Laor et al 1994). Such AGN have FWHM(H/3) 



in the range ~ 500-1500 km/s, which is similar to, or even 
narrower than the FWHM of model A. 

The idea that the inverse quadratic far wings, which 
are observed in many profiles, are linked to parabolic orbits 
and to r~ 2 emissivity dependence, was discussed by several 
authors. The dense and deep potential well in the AGN of- 
fers a 'natural selection' process which favors such orbits 
for the BSs, since all the tightly bound BSs are destroyed 
by collisions or ti dal disruptions. The analytic work of Pen- 
ston et al (1990) assumed strictly parabolic orbits and a 
cloud line emissivity that falls off as r~ 2 . Our numeric re- 
sults confirm that their analysis also holds for a distribution 
of binding energies around E = 0. However, unlike their 
assumed cloud emissivity, the total line emissivity from a 
single BS is nearly constant in r. If we assume that the BSs 
are a fixed proportion of the normal stellar population, the 
resulti ng line profiles are core dominated with very weak 



profile differences. The idea that the BS fraction decreases 
with distance from the black hole is a natural one since such 
stars are not seen in normal stellar environments and their 
existence must therefore be due to the special conditions in 
the AGN. 

It is tempting to interpret the r~ 2 dependence as indi- 
cating that the local photon or particle flux is connected to 
the creation of BSs. However, taken at face value, /bs oc r~ 2 
means only that the BS number is a function of distance 
from the black hole. It implies nothing about the distance 
dependence of the BS properties. In the simplified picture 
presented in this paper, normal stars approach their peri- 
galacton on a near-parabolic orbit, and due to some yet 
unspecified physical mechanism, a small fraction of them 
expands to the BS state at various distances from the black 
hole. If the bloating process occurs in a short time relative 
to the BLR crossing time (e.g. a collisional merger between 
two stars), fss{r) can be interpreted as the probability for 
this to happen to a given star at a given radius. Alterna- 
tively, the bloating may be a gradual process. In this case, 
the simplified two-component stellar population that was 
used here may be thought of as an approximation, where all 
stars bloated beyond a certain threshold are represented by 
a single type of "effective" BS, and the rest by a solar type 
star. /bsM is then the fraction of stars that are bloated 
beyond this threshold at a given radius. 

Although the modified BS model succeeds in reproduc- 
ing the FWs and wide wings to a fair degree, there re- 
main several problems. The calculated profile wings are still 
weaker than the observed ones. This may be due to the fact 
that the cores of the observed profiles have a blended con- 
tribution from the narrow line region. In addition, the dif- 
ferences between the lines are not as distinct as in the mean 
observed profiles. A reliable estimate of this discrepancy is 
possible here only for the C in] A1909/C IV A1549 FWHM ra- 
tio. The Mg II A2798 line is problematic, both in the model, 
which badly under-estimates its strength, and in the ob- 
served spectra, where it is very difficult to deblend. These 
discrepancies may indicate that /bs falls off even faster than 
r~ 2 . 

The question whether BSs can produce the extremely 
broad profiles that are seen in some objects was not inves- 
tigated here. This is part of the larger question of whether 
the BS model has the flexibility to explain the entire range 
of observed AGN spectra, which will be addressed in a fu- 
ture paper. Here, we limit ourselves to the problem of repro- 
ducing the profiles of a typical AGN, where we define typ- 
ical profiles as those that have the median FW and inverse 
quadratic wings. We also did not address the issues of profile 
asymmetry and profile shifts relative to the systemic one. In 
this context, it is interesting to note that there have been 
phenomenological attempts to decompose the BLR emission 
spectrum into several components, and in particular to as- 
sign the line shifts and profile a s ymmetry to a ver y broad 
component dFrancis et al. 1992|; Wills et al. 199$ |Biotli- ' 



wings. rhis_Mtmitio^ i «mj3C_«jn i ccta 



sumptii in and setting /bs oc r~^" The resulting profiles are 
only approximately inverse quadratic, because the emissivity 
of the individual lines is not constant in r. The BS popula- 
tion of this model is more concentrated towards the central 
high velocity region. As a consequence, the differences in 
the emissivity of the various lines begin to manifest itself in 



erton et al. 1994a 



Brotherton et al. 1994b 



Baldwin et al 



1995). In such a scenario, the BSs can be identified with the 



symmetric, unshifted, intermediate width component. 

The conclusion that the BSs are likely to be tidally- 
limited and optically thick has an interesting consequence. 
For a given BLR covering factor, CV, (usually ~ 0.1) the 
total number of BSs, Nbs, is 
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4CV /M bh \ 2/3 ,„,v 

Nbs = xliAm^) ■ (25) 

This explains why in all our models, Nbs ~ 5xl0 4 . It also 
implies that if Xtidal and Mbs are similar in different AGN, 
Nbs/Cf is a measure of the black hole mass. 

To summarize, our attempts to model the BLR line pro- 
files with BSs lead us to the following conclusions. 

(i) The BSs can reproduce the typical full widths of 
observed profiles if Li on /M bh is low enough (e.g. Mbh = 
2xlO 9 M0 and L- ian — 3.6xl0 44 erg/s) and an external cut- 
off mechanism limits the size of the BLR. 

(ii) Due to the small size of the BLR, the most important 
effect limiting the size of the BSs is tidal disruption by the 
black hole. 

(iii) In order to reconcile the observed similarity of AGN 
line spectra over a large range of continuum luminosities 
with the idea of tidally limited BSs, it is necessary to as- 
sume that the wind properties (e.g. mass loss rate or wind 
velocity) are correlated with Mbh and Li on - 

(iv) The BSs can produce the observed wide, inverse 
quadratic far wings if their fraction in the stellar popu- 
lation falls off roughly as r~ 2 . Although this is a highly 
concentrated distribution, the collisional mass loss is small 
enough to be consistent with the observations and BS life- 
times. However, the calculated profiles tend to have weaker 
wings than the observed ones and the differences between 
the profiles of the different species are somewhat smaller 
than those actually observed. 
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APPENDIX: L ion /M bh SCALING OF THE WIND 
MODEL 

In paper I we demonstrated that BSs with dense and slow 
winds reproduce the correct line ratios for one specific AGN 
model. The wind boundary of such BSs is either tidally- 
or mass-limited. Unlike the Comptonization boundary, Tic, 
which establishes itself at a constant U irrespective of Li on 
or Mbh, A r (i? t idai) depends on Mbh, but not on Li on and 
therefore U(Rtidai) depends on both Li on and Mbh- Simi- 
larly, Af(7i m ass) depends on the wind properties, but not on 
Li on and therefore U(Rma.ss) depends on L lon . As a result, 
the emission line spectrum of tidally- or mass-limited BSs 
with given wind parameters depends on the specific choice of 
Li on and Mbh- In paper I we assumed for simplicity that all 
the wind properties, except the wind size, are independent 
of the external environment. It is obvious that this assump- 
tion is inconsistent with tidally- limited BS, since AGN of 
very different luminosities are observed to have very similar 
line ratios. This leads us to hypothesize the existence of some 
unspecified physical mechanism that affects the properties 
of the BSs in a way that makes the emission line spectrum 
insensitive to changes in Li on and Mbh- The results of sec- 
tion ^| indicate that in order to reproduce the observed line 
widths, the AGN Li on /Mbh ratio should be smaller than that 
assumed in paper I. In this case the line emission originates 
closer to the black hole and most of it is from tidally-limited 
BSs. We will therefore concentrate on the scaling properties 
of tidally-limited wind models. 

Let Lion' = ALion and M bb = fiM bb describe the trans- 
formation from an AGN model with Li on and Mbh to one 
with Lion' and Mbh'- Equations [j], ^ and ^ and the relation 
U oc L ion /r 2 N yield 

U(r) oc R t a M BS {2 - a)/3 X?r d « r - a ^L ion M bb (a - 2)/:i . (Al) 

where we approximate the BS mass, Mbs, by IMq, the max- 
imal mass in the BS (star and wind). It then follows that in 
the transformed AGN model 
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U'(r) = \fi 



(«-2)/3 



U{r). 



(A2) 



U changes with r, and therefore the integrated line spectrum 
reflects a weighted sum over a range of U values, weighted by 
the BS DF. The simplest option for making the integrated 
spectrum constant in A and /i is to require that U'(r) = 
U(r). A change in the black hole mass and luminosity imply 
also a change in the stellar DF and the BLR size (i.e. the 
integration limits) . We are assuming here that these changes 
have only a small effect on the integrated spectrum. This 
is indeed the case in the examples which are presented in 
section ^ (cf Fig. |) . 

In principle, any of the free parameters in equation A.1 
can be v aried to compensate for the scaling factor in equa- 
We will assume for simplicity that 



tion 



A2 



(a-3)/3 



(A3) 



For the free-fall, a = 1/2, wind models which we study here, 
this reduces to 



(A4) 



If all AGN emit at some fixed percentage of the Eddington 
luminosity, the scaling law is further reduced to 



(A5) 



The physical interpretation of these scaling laws is be- 
yond the scope of this work. Nevertheless, we note that it 
so ha ppe ns that the phenomenological scaling laws (equa- 
tions A4 and A5) take on a particularly simple and sug- 
gestive form in the relevant case of tidally-limited and free- 
falling wind flows. While this may be a clue to the physics 
of the winds, it should be noted that a simple interpretation 
of these scaling laws in terms of the ionizing flux or tidal 
forces runs into an inconsistency with our simplified wind 
structure, which has no r dependence. 

Unlike a systematic search in parameter space, this scal- 
ing prescription is not guaranteed to yield the best fit results. 
However, the calculations in paper I indicate that the inte- 
grated line spectrum is a strong and smooth function of the 
wind parameters. Small changes in the best fit values of the 
mass loss rate or the wind velocity considerably worsen the 
fit to the line spectrum. This implies that the scaling solu- 
tion, or a solution very similar to it, is likely to be the best 
solution in the region of parameter space that is considered 
in this work. 
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